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f-{— ' ! Abstract 

c/3 ■ A new type of Interface Acoustic Waves (lAW) is presented, for 



single-crystal orthotropic twins bonded symmetrically along a plane 
containing only one common crystallographic axis. The effective bound- 
ary conditions show that the waves are linearly polarized at the in- 

'^ . terface, either transversally or longitudinally. Then the secular equa- 

tion is obtained in full analytical form using new relationships for the 
O ' displacement-traction quadrivector at the interface. For Gallium Ar- 

senide and for Silicon, it is found that the lAWs with transverse (resp. 
longitudinal) polarization at the interface are of the Stoneley (resp. 

^ ; leaky) type. 
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Figure 1: Twinned crystal 



1 Introduction 



The possibilities of existence of Stoneley waves propagating at the interface of 
two welded half-spaces are quite restricted, especially when the semi-infinite 
bodies are made of dissimilar crystals. For example, Owen (1964) famously 
studied 900 combinations of isotropic materials and found only 31 pairs sup- 
porting Stoneley waves. Similar conclusions are also reached when the ma- 
terials are anisotropic. In general, the existence of an interfacial Stoneley 
wave is highly sensitive to the differences in material parameters for each 
medium, in particular to the difference in shear wave speeds (Chadwick & 
Currie, 1974). 

The situation is somewhat more favorable when the half-spaces are made 
of misoriented but identical crystals. For instance, Stoneley waves were found 
to exist for any angle of misorientation and for propagation along the twist- 
angle bisectrix in the case of some hypothetical crystal by Lim & Musgrave 
(1970), Copper by Tholen (1984), Gallium Arsenide by Barnett et al. (1985), 
or Silicon by Mozhaev et al. (1998); they also exist for 180 degrees domain 
boundaries of Barium Titanate or Quartz (Mozhaev & Weihnacht, 1996, 
1997), for propagation in any direction in the interface plane. In these arti- 
cles, the interface is normal to a crystallographic axis which is itself common 
to each half-space. 

Now consider the interface described in Figure 1, where the upper and 
lower half-spaces are made of the same crystal with at least rhombic symme- 
try and with misoriented crystallographic axes X and Y (represented on the 





Figure 2: Cutting (a), rotating (b), and bonding (c) of a rhombic crystal 



Figure making an angle ±9 with the interface and its normal) and common 
crystallographic axis Z (normal to the plane of the Figure). This geome- 
try of chevron or "herringbone" pattern might be encountered in crystals 
subjected to a gliding twinning plane, or to the formation of conjugate kink 
bands through compression. Hussain and Ogden (2000) recently proposed 
a theoretical simulation of the plastic deformation associated with the twin- 
ning of crystals, with a possible exploitation in non-destructive testing of 
materials. Another way of producing such twins would be the following. 
Consider an infinite crystal with orthotropic or higher symmetry; cut it in 
two halves so that the plane interface contains one crystallographic axis (xs 
say) and makes an angle 6 with another crystallographic axis {xi say), see 
Figure 2(a); then rotate one half-space by 180 degrees about the normal to 
the interface, see Figure 2(b); finally rebond the half-spaces, see Figure 2(c). 
A thorough review article by Gosele and Tong (1998) presents several exper- 
imental procedures of wafer bonding (also known as direct bonding or fusion 
bonding or "gluing without glue") and its numerous applications, not only 
for Silicon-based sensors and actuators but also in microsystem technologies, 
nonlinear optics, light-emitting diodes, etc. using Gallium Arsenide, Quartz, 
or Sapphire. 

When the boundary conditions take the symmetries of the problem into 
consideration then the interface wave is found to be linearly polarized at the 
interface, either longitudinally or transversally. Once these effective bound- 
ary conditions are established in Section 2 (and Appendix), two correspond- 
ing secular equations are derived explicitly in Section 3 as cubics in the 
squared wave speed. The results are illustrated graphically for Gallium Ar- 
senide and for Silicon. For these materials, it turns out that one of the secular 
equations corresponds to a non-physical supersonic leaky interface wave, and 
that the other corresponds to a subsonic interface acoustic wave having the 



characteristics of a Stoneley wave: decay with increasing distance from the 
interface, speed larger than that of the corresponding surface Rayleigh wave. 

2 Effective boundary conditions 

Consider Stoneley waves traveling with speed v and wave number k at the 
interface of a bimaterial made of two perfectly bonded orthotropic media. 
Both half-spaces are made of the same crystal (mass density p, non-zero 
reduced comphances s'n, ^22, si2 7 •§44, S55, Sgg). However, the normal to the 
interface 0x2 and the direction of propagation Oxi are inclined at an angle 
6 to the crystallographic axes Oy and Ox of the lower (x2 > 0) half-space 
and at an angle —9 for the upper (2:2 < 0) half-space. 

In the lower half-space, the strain-stress relation is Ej^ = Si^aik where 
the reduced compliances Sik are given by (see for instance Ting (2000) or 
Destrade (2003)), 
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(2.1) 



When the mechanical displacement u{xi,X2,X3,t) is modelled as a linear 
combination of the partial modes g^'^K^i+p^'2--"i)\j (say), then the in-plane 
strain decouples from the anti-plane strain, and so does in-plane stress from 
anti-plane stress (Stroh, 1962); also, p is a root of the characteristic polyno- 
mial , 



u^p - 2u3p + U2P - 2uip + uo = 0, 



(2.2) 



where the (real) coefficients Ui are given in terms of the reduced compliances 
and of X = pv^ by (Destrade, 2003), 
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In the upper half-space, 6 is changed to its opposite so that by fl2.ip . sn, 
■S225 ■S12, See remain unchanged whilst Sie and S2e change signs. Consequently, 
by (12. 3p the characteristic polynomial in the upper half-space is 

uj^p'^ + 2w3p^ + uj2P^ + 2uoip + Wo = 0. (2.4) 

It follows that if p is a root of the characteristic polynomial (12.21) for the 
lower half-space, then —p is a root of the characteristic polynomial (12.41) for 
the upper half-space. 

As a consequence of these properties (change of sign across the interface 
for sie, S2e, and for the p's), the following effective boundary conditions apply 
at X2 = (see Appendix), 

Ml = 0-22 = 0, or U2 = 012 = 0. (2.5) 

Following Mozhaev et al. (1998), interface acoustic waves satisfying the first 
condition (I2.5l) i are denoted lAWl, and those satisfying the second condition 
(12.51) 9 are denoted IAW2. The polarization of the waves at the interface is lin- 
ear: transverse for lAWl and longitudinal for IAW2. Note that Mozhaev and 
collaborators obtained similar effective boundary conditions for other types 
of twin boundaries (180° ferroelectric domain boundary in tetragonal single- 
crystals with crystallographic propagation direction (Mozhaev & Weihnacht, 
1995, 1996); Dauphine twins in quartz with boundary coincident with the 
YZ mirror plane and propagation along the Y direction (Mozhaev & Weih- 
nacht, 1997); twisted single-crystal cubic wafers with basal plane boundary 
and propagation along the twist-angle bisectrix (Mozhaev & Tokmakova, 
1994; Mozhaev et al., 1998). 

3 Equations of motion and explicit secular 
equations 

The equations of motion in the lower half-space are written as a first-order 
homogeneous differential system for the in-plane displacement-traction vec- 
tor, 

${kX2) = [Ui{kx2),U2ikx2),ti2ikx2),t22ikx2)r, (3.1) 

where the Ui and tj2 are defined by 

Ui{xi,X2,X3,t) = Ui{kx2)e'^^''^-''^\ ai2ixi,X2,X3,t) = ikti2{kx2)e'^^''^-''^\ 

(3.2) 



Explicitly, the system is given by (Destrade, 2001) 
^' = zN^, N -- 
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where X = pv"^ and (Ting, 2002) 
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Note that in passing from the lower half-space to the upper half-space, only 
re and n2e change signs. 

Ting (2003) proved that for any positive or negative integer n, the matrix 
N" has the structure. 
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with K(") = K(")t, N^") = N2(")^, (3.5) 



see also Currie (1979). Hence the following matrix M^") is symmetric, 
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Now premultiply both sides of (I3.3p i by — -j^ IN" ^ and add the complex 
conjugate quantity to obtain ^ • M*^"^^ + ^ ■ M*^")^' = 0. Assuming that the 
lAWs vanish at great distance from the interface (^(00) = 0), integration 
between 0:2 = and 0:2 = 00 yields. 



£(0) • M(")^(0) = ^(0) • IN"^(0) = 0. 



(3.7) 



This relationship is valid for any type of anisotropy. Because of the Cayley- 
Hamilton theorem, it yields five linearly independent equations for a 6 x 6 
matrix N, and three equations for a 4 x 4 matrix N, as here. Note that for 
surface wave boundary conditions, the tractions are zero at 0:2 = 0, and the 
relationship reduces to (Taziev, 1989): u(0) • K(")u(0) = 0. At n = 1, 2, -1, 
the following expressions are found for M'^"\ 
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(3.10) 
with A = det N and 



a = (1 - nQQX){r] - X) - r^X, b = -r2 + X{r2nQG - rgnse) 
c = -n22 + X{n22nQQ - nl^^), 

d = -X[rl + n22{r] - X)], e = -r2 - n22(?7 - X), 

/ = 2r2rQn26 - rjuee - r^n22 - (?7 - X){n22nQG - uIq). 



(3.11) 



Corresponding to the first effective boundary condition (l2.5P i is the in- 
terface acoustic mode lAWl, for which the quadrivector ^(0) is in the form, 



^(0) = f/2(0)[0,l,a,0]T, 
(say), so that equations 03. 7p read 

Mg) + m£) (a + a) + M^^^aa = 0. 



(3.12) 



(3.13) 



Choosing in turn n = —1, 1,2, an homogeneous hnear system of equations 
follows, 

ah c 1 

X -1 rim a + a =0. (3.14) 

.0 rQ + n2QX -2(ri26 + renee)] [ aa 

Similarly, corresponding to the second effective boundary condition 02.5P 9 
is the interface acoustic mode IAW2, for which the quadrivector ^(0) is in 
the form, 

C(0) = [/i(0)[l,0,0,a]^, (3.15) 

(say), so that equations (13. 7p read 
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Choosing in turn n = —1,1,2, an homogeneous hnear system of equations 
follows. 
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The secular equations for lAWl and IAW2 correspond to the vanishing of 
the determinant of the 3x3 matrices given in f l3.14p and (I3.17P , respectively. 
They are both cubics in X = pv^ . For purposes of comparison, recall that a 
homogeneous bulk wave propagates in the Ox^ direction at a speed Va such 
that X = pvl is root of the quadratic. 



A = det N = 0, 



(3.18) 



and that a surface Rayleigh wave propagates over either half-space along 
Oxi at a speed vr such that X = pv\ is root of the quartic (Currie, 1979; 
Destrade, 2001; Ting, 2002), 
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Figure 3 shows the influence of the twinning angle 6 upon the wave speeds, 
for (a) Gallium Arsenide and (b) Silicon. The respective elastic stiffnesses 
(10^° N/m^) and mass densities (kg/m^) are (Royer & Dieulesaint, 1996): 
cii = 11.88, Ci2 = 5.38, C66 = 5.94, p = 5307, and cn = 16.56, Ci2 = 6.39, 
C66 = 7.95, p = 2329. Similar comments apply to both crystals. The secular 
equation (I3.17P corresponds to a non-physical leaky wave IAW2 with a speed 
(upper full curve) which is always supersonic (always above the curve for 
the speed given by (I3.18P of a bulk shear wave (upper dotted curve)). The 
secular equation (I3.14p corresponds to a Stoneley wave lAWl with a speed 
(lower full curve) which is always larger than the speed of a Rayleigh wave 
propagating in either half-space (always above the curve for the speed given 
by (I3.19P (lower dotted curve)). 



Appendix 



Here the effective boundary conditions (12. 5p for twinned crystals are derived. 
Solutions to the equations of motion (13. 3p are in the form ^{kx2) = 
^^Qtpkx2^ By substitution, ^o is found from the adjoint matrix to N — pi 
as 

C=[ai,bi,Ci,dif, (A.l) 
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Figure 3: Interface wave speeds in twinned crystals: leaky wave (upper solid) 
and Stoneley wave (lower solid); also included: quasi bulk shear wave (upper 
dashed) and Rayleigh wave (lower dashed). 



with 



Qi = [(1 + r2)X - rjlpi + r2rQX + n2fi{ri - X)X, 
k = -Xp] - 2reXpi - [1 + r^ + n(,,{r] - X)]X + r], 
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2r6p- + [r2 - rg - 71(^^(7] - X)\pi + n2(i{r] - X) + r2r^. (A.2) 



Here p is a root of the characteristic polynomial det (N — pi) = 0, which 
is given explicitly by (12. 2p . For the lower half-space, only the roots pi, p2 
with positive imaginary parts are kept in order to ensure decay away from 
the interface 2:2 = 0. Thus the solution for X2 > is 



^{kx2) = Pieoe'"'"''' + /32^o'e*' 



i2^ip2kx2 



(A.3) 



for some constants /3i and /32- 

In the upper half-space X2 < 0, the p's are changed to their opposite 
so that the roots of the characteristic polynomial with negative imaginary 
parts are —pi, — p2- Moreover, rg and ^26 also change signs, so that the 
eigenvectors ^0 are in the form. 



Sq [ '^i; "i) Cj) (^i 



(A.4) 



and the solution for X2 < is 



(A.5) 



for some constants j3i and /32. 

It follows that the continuity of the displacements and of the tractions 
across the interface 0:2 = reads 



Piil + P2eo = Ml + Ml 

This system of equations is re- written as 



where 
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The determinant of the 4x4 matrix in flA.7p is 4 det A det B. It is zero 
when either (a) det A = 0, and then B/3 — B/3 = yields /3 = /3; or (b) 
det B = 0, and then A/3+ A^ = yields ^ = -j3. In case (a), A/3+A^ = 
means A/3 = , that is 

«i(0) = t22(0) = 0; (A.9) 

in case (b), B/3 — B/3 = means B/3 = , that is 



«2(0) = ti2(0) = 0. 



(AlO) 
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List of Figures. 

Figure 1: Twinned crystal. 

Figure 2: Cutting, rotating, and bonding of a rhombic crystal. 

Figure 3: Interface wave speeds in twinned crystals: leaky wave 
(upper solid) and Stoneley wave (lower solid); also included: quasi 
bulk shear wave (upper dashed) and Rayleigh wave (lower dashed). 

Figure 3(a): Gallium Arsenide. 
Legend on graduated horizontal axis: "propagation direction/crystallographic 
axis angle [deg]." 
Legend on graduated vertical axis: "speed [m/s]." 

Figure 3(b): Silicon. 
Legend on graduated horizontal axis: "propagation direction/crystallographic 
axis angle [deg]." 
Legend on graduated vertical axis: "speed [m/s].". 
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